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PREFACE 


This  report  is  the  first  of  the  two-part  final  documentation  of 
Pacific-Sierra  Research  Corporation's  analysis  of  the  large  urban  fire 
environment.  This  part  develops  the  theory  underlying  the  analysis. 
Parameter  studies  and  simulations  of  model  urban  areas  are  presented 
in  the  second  part.  All  work  was  performed  for  the  Federal  Emergency 
Management  Agency  under  contract  EMW-C-0747;  the  technical  monitor  was 
Dr.  David  Bensen. 

The  authors  gratefully  acknowledge  the  many  valuable  suggestions 
of  Dr.  Harold  L.  Brode. 


SUMMARY 


This  report  analyzes  the  velocity,  temperature,  and  pressure  dis¬ 
tributions  of  a  large  area  fire.  The  region  of  interest  is  the  burning 
area.  The  analysis  considers  the  effects  of  massive  heating  by  combus¬ 
tion;  the  results  presented  are  applicable  to  the  lower  surface  layer. 
Predictions  of  the  basic  street-level  environment  can  be  made  and  above¬ 
shelter  air  temperatures  can  be  defined.  Coupled  with  an  analysis  of 
debris  distributions  and  residual  nuclear  radiation  (fallout  or  induced 
activity)  in  an  attacked  city,  the  velocity  and  temperature  predictions 
will  aid  in  defining  the  environment  facing  rescue  forces. 

This  analysis  differs  considerably  from  previous  treatments.  No 
attempt  is  made  to  extrapolate  results  valid  for  a  free-convection 
column.  The  theory  presented  directly  relates  the  heat  release  of  com¬ 
bustion  and  si2e  of  the  burning  area  to  calculated  velocities  and 
thermodynamic  properties.  The  induced  fire  winds  depend  on  the  size 
of  the  fire  and  the  combustion  rate  and  are,  to  leading  order,  inde¬ 
pendent  of  the  properties  of  the  convective  plume. 

The  results  are  restricted  to  large  fires  (small  fires  generate  a 
different  balance  of  forces  than  those  considered  here).  The  induced 
fire  winds  increase  in  intensity  with  heat  release  or  fire  size.  Com¬ 
pressibility  effects  limit  the  maximum  velocities  attained;  however, 
a  continuous  spectrum  of  fire-wind  velocities  is  predicted,  ranging 
from  very  low  levels  (corresponding  to  the  experimental  Flambeau  fires) 
to  high-velocity  winds  such  as  generated  by  the  1943  Hamburg  firestorm. 

It  is  clear  that  the  fire  is  the  forcing  agent,  and  external  in¬ 
fluences  are  not  needed  to  produce  the  devastating  velocities  encountered 
in  the  firestorms  of  World  War  II.  If  the  fire  persists  long,  swirl 
may  slightly  affect  the  high-velocity  winds  predicted  here. 

Sample  calculations  are  presented  for  the  largest  Flambeau  fires 
and  for  the  Hamburg  firestorm.  The  results  agree  with  the  documented 
observations  and  provide  insight  into  the  detailed  hydrothermodynamic 
interactions  occurring  in  a  large  area  fire. 


-vi- 


A  nuclear  weapon  explosion  can  ignite  fires  over  much  larger  areas 
than  those  covered  by  World  War  II  firestorms.  The  theory  developed  in 
this  report  has  been  applied  to  large  fires  such  as  may  occur  in  modern 
urban  areas  subject  to  megaton-yield  explosions.  Those  results  are 
presented  in  Part  II  of  this  report. 
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SYMBOLS 


A,  B  s  dimensionless  constants 

Cp  “  specific  heat  capacity  at  constant  pressure 
e  »  measure  of  numerical  solution  error 
E  =  error  tolerance 

2ft  ”  effective  turbulent  kinematic  viscosities 

3  *  functional  form  for  ground-level  thermal  boundary 
condition 

g  ■  gravitational  acceleration 
H  *  maximum  height  of  flames 
Hc  =  scale  height  of  convection  column 

k 

k  =  reciprocal  of  graybody  radiation  mean  free  path 
k>l ,  k2  =*  effective  turbulent  thermal  conductivities 
K^,  3  dimensionless  heat-diffusion  coefficients 

M2  =*  dimensionless  momentum-diffusion  coefficients 
P  *  perturbation  pressure  (dimensionless) 

§  »=  pressure 

P.  =  ground-level  atmospheric  pressure  in  far  field 

a 

“  leading-order  perturbation  pressure  for  V  <<  1 

,  P  “  successive  values  of  P  in  numerical  solution  iterations 
old,  new 

q,  q^  ■  dimensionless  spatial  distributions  of  volumetric  heat- 
addition  race 

Q  “  volumetric  heat-addition-rate  scale 
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r  ■=  i adial  position  coordinate  (dimensionless1) 

R  -  fire  radius 

•f* 

T  ■  temperature  (dimensionless  ) 

T  *  ground-level  atmospheric  temperature  in  far  field 

*  leading-order  perturbation  temperature  for  v  <<  1 

1* 

u  ■  radial  velocity  (dimensionless  ) 
u^  *=  leading-order  horizontal  velocity  for  v  <<  .1 
U  *  radial  velocity  scale 

X 

v  =  vertical  velocity  (dimensionless') 
v^  leading-order  vertical  velocity  for  v  <<  1 
w  *  generic  jump  condition  variable 

+ 

y  *  vertical  position  coordinate  (dimensionless  ) 

y  *  rescaled  vertical  coordinate 

y  =•  maximum  v-value  considered  in  numerical  solution 
max 

y  =  specific  heat  ratio 
5  a  di  lensionless  constant 

AR  =  difference  between  maximum  and  minimum  convection  column 
radii  in  turning  region 

e  *  reciprocal  of  combustion-zone  aspect  ratio 
V  *  dimensionless  measure  of  volumetric  heat-addition  rate 
p  a  density  (dimensionless') 
p  *  ground-level  atmospheric  density  in  far  fi^ld 

cl 

o  *  dimensionless  constant 
o  ”  Stefan's  constant 
to  =  relaxation  coefficient 


Except  in  Eq.  (3). 
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t  I ■  INTRODUCTION 


Most  previous  analyses  of  free-burning  fires  have  focused  on  the 
physics  of  the  buoyant  plume.  The  theories  are  restricted  to  small 
fires  and  are  applicable  only  well  above  the  flaming  region.  This 
report  considers  a  large  burning  area  and  analyzes  the  strong-inter- 
action  region  in  and  around  the  fire. 

The  buoyant  plume  above  a  small,  free-burning  fire  is  well  described 
by  integral  theories  that  relate  a  sustained  vertical  flow  to  a  weakly 
compressible  gas  state  [Morten,  Taylor,  and  Turner,  1956].  Dynamic 
pressure  forces  are  assumed  negligible,  and  a  weak  entrainment  of  am¬ 
bient  air  is  postulated  at  an  effective  plume  edge.  In  and  around  the 
burning  zone,  however,  such  theories  have  limited  applicability.  The 
combustion  processes  foster  large  changes  in  temperature  and  density. 
Those  changes  create  a  region  of  strong  buoyancy  and  hence  nonnegligible 
pressure,  gradients. 

Further,  analysis  of  measurements  [Cox  and  Chitty,  1980]  shows 
that  entrainment  criteria  fail  to  accurately  describe  the  flow  field 
at  and  just  above  the  fire  periphery.  The  horizontal  velocity  appears 
to  decay  inversely  with  height  above  the  ground,  and  approaches  a 
variation  proportional  to  the  center-line  vertical  velocity  only  well 
above  the  fire.  It  is  the  buoyancy-generated  pressure  forces,  and  not 
diffusive  entrainment,  that  control  the  low-level  induction  of  ambient 
air  into  the  fire.  The  relationship  between  dynamic  pressure  and  in¬ 
duced  fire  winds  is  indicated  in  the  calculations  of  Smith,  Morton,  and 
Leslie  [1975],  though  they  assume  that  density  changes  are  small  (Bous- 
sinesq  approximation)  and  do  not  seek  fine  resolution  near  the  fire. 

'  •  ;  A  model  of  the  strongly  buoyant  flow  generated  in  and  around  a 

large  area  fire  is  developed  in  Sec.  II,  producing  a  more  highly  re¬ 
solved  description  of  the  mechanics  of  fire-wind  generation.  A  "large 
area  fire"  is  defined  here  as  one  of  large  aspect  ratio,  i.e.,  with  much 
*  larger  horizontal  than  vertical  dimensions.  The  analysis  is  thus  appli¬ 

cable  to  fires  ranging  in  size  from  the  muloiple-fuel-bed  Flambeau 
experiments  [Countryman,  1969]  to  city-scale  fires  arbitrarily  larger 
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than  those  occurring  in  World  War  II  [Bond,  1951;  Irving,  1963;  Middle- 
brook,  1981]. 

Attention  Is  focused  on  the  hydrothermodynamic  processes  in  and 
around  the  combustion  zone,  and  an  uncoupled  boundary-value  problem 
describing  the  nonlinear  relationships  between  heat  addition,  strong 
buoyancy,  and  pressure  gradients  is  developed.  The  induced  fire  winds 
are  related  to  the  characteristics  of  the  combustion  zone  through  jump 
conditions  applicable  at  the  fire  periphery.  The  specification  of 
local  boundary  conditions  enables  one  to  predict  fire  winds  without 
extensive  numerical  computations  involving  regions  far  from  the  combus¬ 
tion  zone  [Smith,  Morton,  and  Leslie,  1975],  and  a  much  more  detailed 
view  of  the  burning-zone  physics  is  made  possible. 

In  one  parameter  limit,  the  model  equations  can  be  solved  in  closed 
form.  The  resulting  solution  concisely  describes  the  basic  interchanges 
of  energy  and  momentum  as  well  as  the  role  of  pressure  gradients  in  fire- 
wind  generation.  Volume  heating  due  to  the  combustion  processes  raises 
the  gas  temperature,  lowers  its  density,  and  creates  buoyancy.  The 
pressure  in  and  around  the  fire  zone  is  then  hydrostatically  lowered, 
inducing  the  fire  winds.  That  inflow  is  kinematically  turned  upward, 
forming  the  initial  part  of  a  convection  column  (plume)  over  the  fire. 

This  analysis  predicts  chat  a  nominal  scale  for  the  fire-wind 
inrush  velocity  varies  linearly  with  fire  size  (width)  and  intensity 
(mean  rate  of  heat  release) .  For  the  multiple-fuel-bed  Flambeau  experi¬ 
ments  and  for  the  World  War  II  Hamburg  firestorm,  velocity  scales  of 
about  9  and  18  m/sec,  respectively,  are  predicted.  Those  values  agree 
well  with  reported  data  [Countryman,  1969;  Middlebrook,  1981].  Simula¬ 
tions  of  the  near-fire  velocity  and  thermodynamic  fields  induced  by 
those  events  have  been  performed,  and  are  presented  in  Sec.  IV.  Again, 
computed  results  duplicate  observed  flow  patterns  and  are  consistent 
with  reported  data. 
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II.  ANALYSIS 


The  basic  features  and  geometry  of  a  large  area  fire  system  are 
illustrated  in  Fig.  1.  The  fuel  bed  may  be  either  continuous  or  com¬ 
posed  of  discrete  packets.  However,  it  extends  far  enough  that 
turbulence  breaks  the  flame  envelope,  tending  to  keep  flame  heights 
uniform  [Thomas,  1963].  Ambient  air  is  drawn  into  the  fire  by  a  com¬ 
bination  of  physical  processes  involving  heat  addition,  buoyancy,  and 
pressure  gradients.  Air  is  heated  in  the  combustion  zone  and  rises 
buoyantly,  carrying  the  products  of  combustion,  to  form  a  convection 
column.  As  cooler  air  replaces  the  rising,  heated  air,  a  high-speed 
inflow  is  produced  in  a  surface  layer  that  includes  the  lowest  part  of 
the  convection  column  as  well  as  the  entire  combustion  zone.  The 
portion  of  chat  layer  that  extends  the  width  of  the  combustion  zone  is 
here  called  the  turning  region. 

The  present  analysis  examines  the  hydrothermodynamics  of  tv>e 
turning  region.  In  doing  so,  it  also  provides  the  initial  conditions — 
vertical  velocity  and  density  difference  from  ambient — required  for 
treatment  of  the  convection  column  flow. 

The  turning-region  flow  is  taken  to  be  quasi-steady  and  axisym- 
metric.  In  Fig.  1,  R  represents  the  fire  radius,  H  the  maximum  flame 
height,  AR  the  difference  between  the  maximum  and  minimum  radii  attained 
by  the  convection  column  in  the  turning  region,  and  Hc  a  scale  height 
for  the  column.  For  large  area  fires,  H  <<  R.  A  characteristic 
combust ion- zone  aspect  ratio  t  ^  is  defined  by 

e  -  |  «  1  .  (1) 

It  has  been  observed  that  the  column  radius  of  a  large  fire  de¬ 
creases  just  slightly  with  height  in  the  first  five  to  fifteen  flame 
heights,  then  increases  slightly  with  altitude  [U.S.  Forest  Service, 
1967].  The  initial  decrease  in  radius  can  be  associated  with  the  strong 
fire-wind  convection,  so  the  height  (''pinch  point")  at  which  the  radius 
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first  begins  co  increase  provides  a  reasonable  measure  of  Che  effective 
height  of  the  turning  region.  In  the  model,  it  is  thus  assumed  that 
that  height  is  0(H)  and 


AR  <<  R  .  (2) 

The  net  effect  of  the  combustion  processes  is  modeled  as  a  spatially 
dependent  rate  of  heat  addition.  The  basic  flow  and  thermodynamics  of 
the  turning  region  are  thus  decoupled  from  the  specifics  of  those  pro¬ 
cesses,  and  the  complexities  of  detailed  combustion  modeling  avoided. 
Changes  in  gas  composition  during  combustion  are  also  assumed  negligible, 
but  the  Boussinesq  approximation  is  not  employed  because  large  changes 
in  temperature  and  density  are  expected.  The  flow  in  the  turning  region 
is  thus  taken  to  be  that  of  an  ideal,  compressible  gas  undergoing  heat¬ 
ing  in  a  finite  zone  at  the  bottom  of  the  region. 

There  is  some  conjecture  that,  after  long  periods  of  time,  ambient 
vorticity  may  be  sufficiently  concentrated  to  engender  a  rotating  or 
swirling  column  [Long,  1967].  The  rotation  would  in  turn  apply  a  biasing 
radial  pressure  gradient  on  the  burning  region  and  thus  affect  the  rate 
of  airflow  into  the  fire.  However,  only  a  few  observations  support  that 
hypothesis  [Carrier,  Fendell,  and  Feldman,  1982].  Therefore,  rotational 
forces  are  neglected  in  our  analysis. 

The  experimental  work  cL  McCaffrey  [1979]  and  the  analyses  of 
Murgai  [1962]  and  Smith  [1967]  suggest  that  radiation  from  the  heated 
gas  and  attendant  smoke  plays  an  important  role  in  reducing  local  air 
temperatures  to  nearly  atmospheric  within  a  few  flame  heights  above 
the  fire.  Radiation  losses  as  well  as  the  effects  of  turbulent  mixing 
are  therefore  included  in  the  model,  though  in  rather  simple,  qualita¬ 
tive  ways.  The  radiation  losses  are  assumed  to  occur  in  volume  and  to 
be  graybody,  or  "transparent"  [Murgai,  1962].  Eddy  dif fusivities  are 
used  to  describe  the  Reynolds  stresses. 

At  the  fire  periphery  (radius  »  R) ,  steep  gradients  in  temperature, 
pressure  [Smith,  Morton,  and  Leslie,  1975],  and  mixing  coefficients  are 
to  be  expected.  Though  they  occur  over  distances  much  less  than  R,  such 


•'«e»w>n<»t>nWK  assay— a? wi  »a 
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rapld  changes  are  critical  in  determining  the  more  gradual  variations 
over  the  entire  width  of  the  turning  region.  To  account  for  them  in 
the  model,  the  following  assumptions  are  also  made.  Since  AR  <<  R,  the 
radius  of  the  turning  region  is  taken  to  be  R  at  all  altitudes,  and 
rapid  changes  in  quantities  around  the  periphery  are  idealized  as  jumps 
(discontinuities)  at  that  radius.  The  magnitudes  of  such  jumps  are 
governed  by  standard  statements  of  mass,  momentum,  and  energy  conser¬ 
vation.  Outside  the  fire,  the  thermodynamic  state  is  assumed  to  be  that 
of  the  ambient  atmosphere  and  turbulent  stresses  are  taken  to  be  much 
less  than  those  inside. 


SCALINGS  AND  BASIC  EQUATIONS 

We  seek  leading-order  predictions  of  the  mean  horizontal  and  ver¬ 
tical  velocity  u  and  v,  temperature  T,  density  p,  and  pressure  P  in  the 
turning  region.  The  following  nominal  variable  scalings  (denoted  by 
braces)  are  employed: 


(r)  -  R  , 

{y}  *»  H  , 

{u}  -  U  , 

{v}  -  eU  , 

{T}  *  T  , 

3 

{p}  -  Pa  . 

(p>  -  pa , 

(3) 

where  r  and  y  represent  radial  and  vertical  position,  respectively; 

U  is  a  horizontal  velocity  scale  yet  to  be  chosen;  and  T  ,  p  ,  and  P 

3  3  3 

represent  ground-level  atmospheric  temperature,  density,  and  pressure. 
Since  a  subsonic  flow  is  expected,  we  rescale  pressure  using 


A 

f-  “  1  +  <$P  , 


(4) 


and  anticipate  that  6  <<  1.  Subject  to  this  rescaling,  the  leading- 
order  set  of  conservation  and  state  equations  is  then  found  [Small, 
Larson,  and  Brode,  1981]  to  be 


1 


t  i 


O 
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~  (rpu)  +  ~  (rpv)  -  0  , 


/  Jui  .  3u\  3P  ,  u  / 1  3  /  3u\  u  \  „  32u 

\  3r  +  V  3y  j  "  '  3r  +  M1  (  r  3r  (r  Sr  )  “  +  M2  ^ 2 


(“  I?  ■■■ w  U)  ■  B(q(c' 


y)  -  cr(T  -  1) 


) 


•■>(«(■  8)) 


+  K, 


2  3y2 


PT  -  1  , 


I 


(5a) 

(5b) 


(5c) 


(5d) 

(5e) 


where 


Si 

a2 


B 


^2 

PaUH  • 

M2  "  epaUH  ’ 

ekj^ 

k2 

rr  n 

cPPaUH  ’ 

2  cpepaUH 

Airak*  ~ 

■  (it)  • 

(6) 


Here,  Q  is  a  scale  for  the  rate  of  heat  addition  due  to  combustion  and 
q(r,  y)  a  specified  spatial  distribution.  The  and  k.^  are  dimen¬ 
sional  mixing  coefficients,  the  specific  heat  capacity  cp  is  assumed 
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constant,  6  ia  Stefan's  constant,  and  k  is  the  reciprocal  of  the  radia¬ 
tion  mean  free  path  (assumed  constant). 

An  appropriate  value  for  the  horizontal  velocity  scale  U  is  found 
by  balancing  the  terms  for  convective  transport  and  heat  addition  in 
the  energy  equation  90  as  to  properly  represent  the  physics  of  a  flow 
driven  by  combustion  heating.  Accordingly,  we  take  B  -  1  and  find  the 
characteristic  fire-wind  velocity  to  be 


For  the  sample  fires  considered  in  Sec.  IV,  typical  areal  heating  rates 

(QH)  and  burning-zone  aspect  ratios  (e  are  on  the  order  of  5  *  10^ 

2 

cal/m  -sec  and  20,  respectively.  The  characteristic  velocity  scales  are 
thus  of  order  10  m/sec. 

BOUNDARY  CONDITIONS 

The  type  of  boundary  conditions  to  be  used  with  Eqs.  (5)  depends 
on  the  basic  nature  of  the  energy  and  horizontal  momentum  equations, 
and  hence  on  the  relative  magnitudes  of  the  coefficients  and  K^, 
i  *  1,  2.  A  general  model  of  turning-region  turbulence  might  assume 
all  such  coefficients  to  be  of  numerical  order  one.  The  energy  and 
horizontal  momentum  equations  would  then  be  elliptic.  Special  cases 
arise,  however,  when  M^,  <<  M,,  ~  0(1)  or  M^,  <<  M^,  0(1). 

Equations  (5b)  and  (5d)  are  then  nearly  parabolic,  and  it  is  more  appro¬ 
priate  to  consider  a  parabolic  analysis  than  a  fully  elliptic  one.  In 
this  subsection,  we  discuss  the  boundary  conditions  required  for  each 
type  of  analysis. 

We  consider  first  the  elliptic  case,  where  all  mixing  coefficients 
are  of  order  one.  The  flow  must  satisfy  a  no-slip  condition  at  the 
ground,  and  a  temperature  condition  must  also  be  enforced.  That  is, 


u  *  v  ■  0 


3(T>  !?)• 0 


on  y  *  0  , 


(8) 


» 
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where  3-0  prescribes  either  a  heat  loss  to  the  ground  or  a  temper¬ 
ature  distribution.  Symmetry  at  the  axis  requires 

u  "  3r  "  0  on  r  *  0  ,  (9) 

For  the  flow  in  the  upper  part  of  the  turning  region  to  be  of  the 
nearly  vertical,  weakly  buoyant  type  characteristic  of  the  convection 
column,  it  is  further  necessary  that,  to  leading  order, 

u  -►  0  ,  T  -►  1  ,  P  +  Ay  -*■  0  as  y  .  (10) 

The  conditions  in  Eqs.  (10)  are  formally  derived  [Small,  Larson,  and 
Brode,  1981]  from  a  matching  of  asymptotic  expansions  that  describe, 
as  e  -*•  0,  che  separate  turning-region  and  convection-column  flows. 

Since  H  is  much  less  than  the  atmospheric  scale  height,  we  take  the 
ambient  temperature  and  pressure  outside  the  turning  region  to  be 
Ta(l  +  o(l))  and  Pa(l  +  6Ay(l  +  o(l))),  respectively  [see  Eqs.  (3) 
and  (4)]. 

Finally,  in  order  to  conserve  mass,  momentum,  and  energy  across 
the  fire  perimeter  (r  ■  1),  the  magnitudes  of  effective  jumps  in  physi¬ 
cal  quantities  at  r  *  1  must  satisfy  certain  conditions.  Those  jump 
conditions  are  easily  found  by  writing  Eqs.  (5a),  (5b),  and  (5d)  in 
conservation  form  and  then  integrating  them  from  r  *  1  to  1+.  With  the 
jump  (r  =  1  to  1  )  being  denoted  by  [w],  the  resulting  conditions  are 

[pu]  -  0  , 

IPUZ]  -  -  [P]  +  k  If]  , 

■  [h  w]  ■  (11) 


Since  Eq.  (5c)  contains  no  r-derivative ,  the  analogous  jump  con¬ 
dition  is  trivially  satisfied  and  a  consistent  leading-oi.de r  description 
is  maintained. 
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To  leading  order,  with  mixing  coefficients  being  much  greater  in¬ 
side  the  turning  region  than  outside,  and  atmospheric  temperature  being 
T^(l  +o(l)),  we  have 


+  1 

r«  5ui 

/3u\“ 

r„  ail 

v-  /  3T\“ 

-  p  ■  1  » 

LM1  3r. 

“  M1 

\3r )  * 

"  K1  (  37/ 

i  4* 

where  T  ■  T  .  ,  etc.  Thus,  by  integrating  Eq.  (5c)  along  r  *  1  , 

1  r«*l 


P  *  -Av  . 


(13) 


Using  Eqs.  (12)  and  (13),  we  can  rewrite  the  jump  conditions  in 
Eqs.  (11)  as 


u  -  pu  , 

W  du  /n  i  a  \  -  2  ,  +.2 

M1  Jr  ~  (p  +  Ay)  *  pu  "  (u  )  > 

ST  *4* 

K1  3r  u  **  u  »  (14) 

where  nonsuperscripted  variables  represent  values  at  r  «  1  .  Substi¬ 
tuting  Che  first  equality  here  into  the  second  and  the  third,  then 
using  Eq.  (5e) ,  we  find  the  following  boundary  conditions  to  be  appli¬ 
cable  at  the  fire  periphery: 

M1  !?  *  (?)”  <T  -  1)  +  (P  +  Ay)  ,  Kl  g  -  f  (T  -  1) 


on  r  «  1  .  (15) 

The  primarily  elliptic  boundary-value  problem  posed  by  Eqs.  (5), 
(8),  (9),  (10),  and  (15)  provides  a  model  for  the  hydrothermodynamics 
of  the  turning  region  in  the  case  where  all  and  are  of  numerical 
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order  one.  That  model  thus  describes  flows  for  which  the  effects  of 
turbulent  mixing  are  dynamically  important  in  both  horizontal  and  ver¬ 
tical  directions. 

We  now  consider  the  nearly  parabolic  case  where  M^,  K^, 

which  is  characteristic  of  flows  having  a  vertically  dominant  turbu¬ 
lent  structure.  In  a  fully  parabolic  analysis  with  M  ■  0, 

the  boundary  conditions  in  Eqs.  (8)  and  (10)  must  be  retained,  but  the 
symmetry  conditions  in  Eq.  (9)  cannot  be  imposed.  In  general,  those 
symmetry  conditions  can  be  satisfied  only  through  the  development  of  an 
axial  boundary  layer  in  which  horizontal  mixing  again  plays  a  leading- 
order  role.  Additionally,  the  two  jump  conditions  in  Eqs.  (15)  degenerate 
for  ■  0  to  the  single  condition  T  *  1.  With  no  peripheral  con¬ 

straint  being  imposed  on  u,  the  parabolic  analysis  cannot  yield  a  unique 
determination  of  the  turning- region  flow. 

The  other  nearly  parabolic  case  occurs  for  flows  having  a  hori¬ 
zontally  dominant  turbulent  structure,  i.e.,  where  M^,  <<:  K^. 

A  developing  column  flow  has  such  a  structure.  In  a  fully  parabolic 
analysis  with  ^  "  0,  all  boundary  conditions  in  Eqs.  (9)  and  (15) 

must  be  retained,  so  as  to  preserve  the  correct  flow  structure  about  the 
symmetry  axis  and  at  the  fire  perimeter.  However,  along  y  ■  0,  neither 
the  no-slip  condition  on  velocity  nor  any  condition  on  temperature  can 
be  enforced.  The  satisfaction  of  those  conditions  [Eqs.  (8)]  requires 
the  development  of  a  thin  surface  layer  in  which  vertical  mixing  plays 
an  important  role.  Further,  the  only  necessary  condition  in  Eqs.  (10) 
is  that  involving  pressure.  As  Eqs.  (5c)  and  (5e)  imply,  p  and  T  -*•  1 
as  y  +  °“  if  P  +  Ay  ■+■  0.  Thus,  from  Eqs.  (5b),  (9),  and  (15),  u  -*•  0 
in  that  limit.  The  proper  set  of  boundary  conditions  for  the  case 
where  M2  -  *  0  is  therefore 

v  -  0  on  y  *  0  ,  (16a) 


9T  n 
u  *  97  -  0 


(16b) 


'  T" 
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P  +  Ay  «  Q 


(16c) 


i  3u  /u\2  , 

1  3r  *  (t)  (T  “  1)  +  (P  +  Ay)  , 


f  (T  -  1) 


on  r  «  1 


(16d) 


~  rzr::,;i„-  r.  ■rr:r~r  ~  - 
s  r;  rr  “err 

h  resulting  boundary-value  p.oble.  is  .  local  one;  Us  so 

r ;  rdlctio° °f  cha  near-eir=  -  *. 

fire  winds  without  additional  analysis  of  the  £„  field.  rhe  ^ 
actions  present  soiutlons  and  results  for  that  parabolic  case. 
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III.  SOLUTIONS 


CLOSED-FORM  SOLUTION  FOR  WEAKLY  HEATED  FLOWS 

Before  considering  the  general  numerical  solution  of  Eqs.  (5)  and 
(16)  for  M2,  Kj  “  it  is  instructive  to  examine  the  restricted  case 
of  a  "weakly  heated"  flow,  i.e.,  q(r,  y)  small.  A  representative 
closed-form  solution  that  summarizes  the  physical  interactions  in 
the  turning-region  flow  is  derived  using  the  following: 


q(r,  y)  -  vq^y)  , 


qx(y) 


for  0  sy 
for  y  >  1 


\)  «  1  . 


(17) 


For  weak  heating,  relatively  small  temperature  changes  and  velocities 
are  expected.  Asymptotic  expansions  of  the  following  type,  valid  as 
v  0,  are  thus  sought  for  the  solution  of  Eqs.  (5)  and  (16): 


T  ~  1  +  VTj.  +  •••  ,  (18a) 

p  ~  1  -  VT^  +  •••  ,  (18b) 

P  ~  -Ay  +  vPx  +  •••  ,  (18c) 

u  ~  vu^  +  • ■ •  ,  (18d) 

v  ~  vv^  +  • » •  .  (18e) 


In  the  V  -*■  0  limit,  the  leading-order  problem  obtained  by  substituting 
Eqs.  (17)  and  (18)  into  Eqs.  (5)  and  (16)  is  then 
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subject  to 


v  »  0 


on  y  ■  0  , 


■  .2i 

1  3r 


on  t  »  0  , 


?L  -v  0 


as  y 


3u! 

*1  TT  "  Pi  - 


Hi 

Sr 


on  r  =  1  , 


(19a) 

(19b) 

(19c) 

<19d) 

(20a) 

(20b) 

(20c) 

(20d) 


The  problem  posed  by  Eqs.  (19)  and  (20)  is  readily  solved.  Equa¬ 
tion  (19a)  involves  ^  alone.  Subject  to  boundary  conditions  (20b)  and 
(20d),  it  must  have  a  solution  that  depends  only  on  y.  Once  that  solu¬ 
tion  is  found,  Eqs.  (19b),  (19c),  and  (19d)  are  decoupled  and  solved  in 
succession  for  u^  and  v^  For  y  >  1,  q^y)  =  0.  The  only  solu¬ 
tion  of  Eqs.  (19)  then  satisfying  Eqs.  (20b),  (20c),  and  (20d)  has  T 
Pl»  ui>  and  Sv^/Oy  all  zero.  That  is, 

T1  =  P1  =  U1  5  0  *  v]_  =  vi(r.  1)  for  y  >  1  ,  (21) 

where  the  v^r,  1)  profile  is  yet  to  be  found. 
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For  0  S  y  S  1,  q  (y)  5  1  and  it  can  be  easily  verified  that  the 
unique  solution  of  the  uncoupled  boundary-value  problem  for  is 

=  l/(4o).  Formally,  then,  the  solution  for  suffers  a  discontin¬ 
uity  at  y  ■  1.  However,  the  discontinuity  may  be  removed  by  using  a 
boundary-layer-type  solution  about  y  *  1.  With  y  stretched  in  this 
layer  by 


(22) 


(23) 


An  appropriate  solution  of  this  system  provides  the  desired  smooth 
transition  in  from  1/  (4c)  as  y  *  -°°,  to  zero  as  v  +».  It  is  fur¬ 
ther  clear  from  Eqs.  (23)  that  and  v^  must  be  independent  of  y  over 

the  boundary  layer.  The  solution  of  Eqs.  (19)  for  0  i  y  S  1  must  there¬ 
fore  satisfy  P^(r,  ^  and  provides  the  v^(r,  1)  profile  to  be  used 

for  y  >  1. 

With  5  l/(4o)  and  P^(r,  1)  *  0,  the  complete  solution  of  Eqs. 
(19)  and  (20)  for  the  fire  zone  is  then  found  to  be 


r  . 
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P1  “  ‘  to  (1  ‘  y)  * 

/Pl\  A  ..  . 

U1  m{\)X  "  “  toM^  (1  '  y)r  ’ 

V1  "  2^(y  "  Y)  for  0  <  y  s  1  .  (24) 


This  solution  summarizes  the  basic  turning-region  physics.  The  com¬ 
bustion  heating  [Eq.  (19a)]  causes  an  increase  in  temperature  and  thus 
a  decrease  in  density  [Eq.  (18b)].  The  mean  pressure  is  then  hydro¬ 
statically  lowered  [Eq.  (19b)],  inducing  the  fire-wind  inflow  [Eq. 

( 20d ) ] .  Finally,  the  inflow  is  kinematically  turned  upward  [Eq.  (19d)], 
forming  the  initial  part  of  the  convection  column. 

NUMERICAL  SOLUTION  PROCEDURE 

In  general,  the  boundary  value  problem  posed  by  Eqs.  (5)  and  (16) 
must  be  solved  by  numerical  computation.  Subject  to  M,,  “  °»  the 

differential  equations  in  set  (5)  are  either  parabolic  or  first-order 
hyperbolic.  Such  equations  are  usually  solved  subject  to  initial  con¬ 
ditions  or  certain  combinations  of  initial  and  boundary  conditions,  but 
not  subject  to  boundary  conditions  on  all  sides  of  the  domain  of  defini¬ 
tion,  as  prescribed  by  Eqs.  (16).  The  boundary  value  problem  is  thus 
unusual,  and  a  novel  algorithm  is  employed  to  solve  it  numerically. 

Our  scheme  involves  iteration  toward  the  correct  form  of  P(r,  0) 
[Small,  Larson,  and  Erode,  1981].  For  y  **  0,  we  have  v  =  0,  and  u,  T, 
and  P  must  satisfy  the  following  ordinary-differential-equation  boundary 
value  problem: 
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(?)fe--E*«i(i£(rfe)-2)  ■ 

(l)g-q(r,  0)  -o«4-l)+K1(i^(r|i))  , 


M 


3u 
1  3r 


(T  -  1)  +  P  , 


f  (T  -  1) 


at  r  *  1  , 


(25) 


First,  a  form  is  chosen  for  P(r,  0),  then  u(r,  0)  and  T(r,  0)  are 
predicted  from  a  standard  finite-difference  solution  of  Eqs.  (25).  A 
nonlinear  Crank-Nicolson  scheme  [Isaacson  and  Keller,  1966]  is  next 
used  to  obtain  a  numerical  solution  of  Eqs.  (5),  (16a),  (16b),  and  (16d) 
for  y  >  0.  After  the  calculation  is  carried  to  some  specified  altitude 
y  (>>H),  the  error  index  e  ■  ]P(r,  y  )  +  Ay  I  is  computed 

[cf.  Eq.  (16c)].  If  e  <  E,  with  E  being  a  specified  error  tolerance, 
the  solution  procedure  is  completed.  If  e  >  E,  a  new  form  for  P(r,  0) 
is  chosen  using  the  following  simple  iteration  scheme: 

P(r’  0)new  "  P(r>  0)old  "  U  (P(r*  W  +  Aymax)  ’  (26) 

where  u  is  a  specified  relaxation  coefficient.  The  prediction  steps 
are  then  repeated  until  the  process  converges  with  e  <  E. 

In  practice,  the  scheme  has  proved  to  be  stable  and  convergent — 
subject,  of  course,  to  reasonable  initial  values  for  P(r,  0).  For 
q(r,  y)  as  in  Eqs.  (17),  the  resulting  numerical  solution  reproduces 
the  closed-form  weak-heating  solution  in  Eqs.  (18),  (21),  and  (24). 


IV.  RESULTS 


FLAMBEAU  FIRES 

The  aultiple-fuel-bed  Flambeau  fires  [Countryman,  1969;  Palmer, 
1981]  were  large,  controlled  burns  conducted  to  simulate  the  character 
istics  of  a  fire  that  has  simultaneously  engulfed  nearly  all  parts  of 
an  urban  area.  The  fuel  beds  consisted  of  several  hundred  separated 
piles  of  mixed  juniper  and  pinyon  pine  logs  and  cuttings.  Each  experi 
ment  was  instrumented  for  temperature,  velocity,  and  weight  loss. 
Additionally,  the  fires  were  filmed  to  provide  a  permanent  visual 
record. 

The  largest  Flambeau  fires  covered  areas  of  about  20  ha,  had 

flames  rising  up  to  20  m,  and  produced  heat  release  rates  on  the  order 
4  5  2 

of  10  to  10  cal/m  -sec  [Countryman,  1969;  Palmer,  1981].  Maximum 
heat  release  is  chought  to  have  occurred  in  and  around  the  fuel  zone 

(height  '-4  m),  the  heating  decreasing  with  altitude  above  that  zone. 

j- 

In  simulating  the  largest  Flambeau  fires,  we  thus  use 


R  -  250  ra  ,  H  -  20  m  ,  QH  -  5.7  x  104  cal/m“-sec  , 


j  1.6 

for 

0  £  y  £0.25 

q(r,  y)  *  1.6  (1  -  y)j 

for 

0.25  £  y  £  1.0 

(o 

for 

y  S:  1.0 

(27) 


From  Eqs.  (1),  (4),  (6),  (7),  and  (27),  we  then  have  e  =>  0.08,  a  hori¬ 
zontal  fire-wind  velocity  scale  of 


The  specific  value  of  QH  here  is  suggested  by  Lommasson  et  al. 
[1968].  The  specific  form  of  q(r,  y)  is  chosen  to  satisfy  the  normal¬ 
ization 

1 

rq  (r ,  y)  dr  dy  -  1/2  , 


which  holds  if  q  =  1. 


» 


U  ■  8.41  m/sec  , 


6  -  9.04  x  10~H  ,  A  -  2.77  .  (29) 

The  required  data  set  is  completed  by  specifying  a  radiation  mean 

free  path  k*'’1  and  turbulent  mixing  coefficients  M  and  K  .  Use  of  the 

1 

graybody  approximation  [Murgai,  19621  assumes  that  k  H  «  1.  However, 
Smith  [1967]  suggests  that,  for  a  turbulent  flow,  the  effective  absorp¬ 
tion  coefficient  may  be  30  to  80  times  larger  than  the  laminar  absorption 

★ 

coefficient,  implying  that  k  H  could  be  of  order  one  or  larger.  Since 

measurements  supporting  estimates  of  that  parameter  are  not  available 

* 

for  the  fires  considered,  we  select  the  nominal  value  k  H  =  1.  Thus, 
from  Eqs.  (27)  and  (6), 

k  »  20  n  and  a  =  0.022  .  (30) 

The  one-parameter  eddy-viscosity  model  of  turbulence  effects  represents 
a  gross  simplification  of  the  actual  field  turbulence.  Because  of  a 
lack  of  reliable  correlations  or  other  relevant  data,  more  detailed 
modeling  does  not  seem  warranted  at  this  time.  For  that  reason,  we 
have  performed  calculations  for  several  sets  of  eddy  coefficients.  The 
sample  results  we  now  present  were  developed  for 

M  -  Kx  -  4.0  ,  (31) 

which  corresponds  to  S./p  U  »  1000  m. 

-L  3 

The  results  of  the  Flambeau  simulation  are  illustrated  by  the 
vector,  contour,  and  profile  plots  in  Figs.  2  through  8.  The  basic 
flow  pattern  predicted  by  the  model  analysis  is  diagrammed  in  Fig.  2. 

The  temperature  rise  and  subsequent  pressure  drop  in  and  around  the 
fire  due  to  combustion  heating  are  shown  in  Figs.  3  and  4,  respectively. 
The  temperature  attains  a  maximum  mean  of  637°K  at  the  center  of  the 
combustion  2one.  However,  it  then  decreases  rapidly  with  increasing 
height,  attaining  a  state  of  weak  buoyancy  above  an  altitude  of  100  ra. 
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Fig.  2--Velocity  field  for  multiple-fuel-bed  Flambeau  fire 
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Pressure  is  almost  independent  of  radial  range  except i  of  course,  at 
the  turning-region  periphery.  Such  behavior  is  also  characteristic  of 
the  more  global  computational  results  presented  by  Smith,  Morton,  and 
Leslie  [1975].  The  maximum  pressure  drop  at  the  periphery  occurs  at 
ground  level,  and  is  0.060  psi. 

Figure  5  illustrates  the  basic  nature  cf  the  horizontal  inrush 
induced  by  the  overall  pressure  drop.  The  maximum  fire-wind  inflow 
velocity  at  the  periphery  also  occurs  at  ground  level,  and  is  9.29  m/sec. 
The  resultant  upward  flow  consistent  with  the  induced  inflow  is  described 
in  Fig.  6.  Vertical  velocity  profiles  are  all  nearly  "top-hat,"  with 
v  «  4  m/sec  in  the  upper  part  of  the  turning  region. 

Figure  7  plots  the  center-line  variation  in  temperature.  In  the 
heating  zone,  the  temperature  reaches  a  maximum  and  then  decreases 
owing  to  radiative  cooling  of  the  hot  gases.  As  y  -*•  °»,  the  temperature 
asymptotes  to  that  of  the  ambient  atmosphere.  The  center-line  velocity 
variation  in  Fig.  8  shows  a  rapid,  nearly  linear  increase  in  the  turning 
region,  then  a  slower  increase  with  altitude.  For  y  >,  7,  the  velocity 
approaches  a  constant  value  consistent  with  the  initial  convection 
column  velocity.  This  behavior  is  similar  to  that  observed  for  small 
fires  [McCaffrey,  1979]. 

In  general,  quantitative  comparisons  of  model  predictions  with 
Flambeau  measurements  are  hindered  by  measurement  scatter  and  the  proba¬ 
bility  of  error.  However,  a  number  of  qualitative  comparisons  can 
realistically  be  made;  in  those  cases,  the  predictions  match  the  ob¬ 
served  physical  phenomena.  For  example,  the  basic  flow  pattern  in  Fig. 

2  is  consistent  with  filmed  observations  of  the  largest  Flambeau  fires 
[U.S.  Forest  Service,  1967]:  the  horizontal  inflow  is  turned  strongly 
upward,  with  a  nearly  vertical  convection  column  flow  beginning  to 
form  at  several  flame  heights  above  the  fire.  Additionally,  predicted 
inflow  velocities  at  the  periphery  and  predicted  temperatures  match 
well  with  (means  of  the)  values  reported  for  those  fires  [Countryman, 
1969].  The  peripheral  fire-wind  profile  in  Fig.  5  is  also  consistent 
with  the  filmed  observations:  below  the  pinch-point  level  (see  Fig.  1), 
maximal  inward  smoke  advection  occurs  at  ground  level,  but  decreases 
with  increasing  altitude. 


o  0.1  0.2  0.4  0.6  0.8  1.0  ’ 

Radius  {  m  x  250  ) 

Fig.  5--Horizontal  velocity  profiles  for  multi  pie-fuel -bed 

Flambeau  fire 
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Fig,  6--Vertical  velocity  profiles  for  nultiple-fuel-bed 

Flambeau  fire 
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|  HAMBURG  FIRESTORM 

The  firestorm  that  devastated  the  Hamoerbrook  and  Borgfelde  dis¬ 
tricts  of  Hamburg  in  July  1943  [Middlebrook,  1981]  was  a  result  of  un¬ 
usually  concentrated  incendiary  bombing.  Those  districts  were  largely 
|  residential,  with  densely  packed,  three-  to  six-story  wood  and  masonry 

buildings  providing  the  primary  fuel  bed.  More  technical  data  are 
available  for  this  event  than  for  other  World  War  II  firestorms;  never¬ 
theless,  the  data  are  limited. 

;  In  simulating  the  peak-period  firestorm,  we  use 

R  ■  1500  m,  H  *  60  m  ,  QH  *  5.7  *  10^  cal/m^-sec  (32) 

[Lommasson  et  al.,  1968;  DCPA,  1973;  Middlebrook,  1981].  Additionally, 
since  the  Hamburg  fuel-zone  height  ~  15  m,  we  take  q(r,  y)  to  be  as  in 
Eqs.  (27).  From  Eqs.  (1),  (4),  (6),  (7),  and  (32),  we  then  have 
e  *>  0-04,  a  fire-wind  velocity  scale  of 

U  *  16.8  m/sec  ,  (33) 

and 

6  -  3.61  x  10"3  ,  A  -  2.08  .  (34) 

*-l 

Consistent  with  our  previous  assumptions,  we  assume  k  »  20  m 
[cf.  Eqs.  (30)]  and  thus  from  Eqs.  (6)  find 

* 

0  =  0.066  .  (35) 

Similarly,  we  performed  model  calculations  for  various  sets  of  mixing 
1  coefficients  and  present  results  here  for 

M1  -  K  *  2.0  .  (36) 

3.  In  general,  the  value  of  the  eddy  viscosity  should  reflect  the  eddy 

size  as  wall  as  the  magnitude  of  the  mean  flow.  Intuitively,  it  is 
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expected  that  the  eddy  size  should  increase  with  the  fire  radius.  The 
change  in  parameter  values  between  Eqs.  (31)  and  (36)  corresponds  to 

a  linear  increase  in  &,/p  with  radius. 

1  a 

The  results  of  the  Hamburg  simulation  are  illustrated  in  Figs.  9 
through  13.  The  temperature  contours  plotted  in  Fig.  10  show  the 
effect  of  the  combustive  heating  occurring  in  the  burning  zone.  At 
the  fire  periphery  (r  ■  1,  y  £  1) ,  the  temperature  jumps  from  the  am¬ 
bient  to  a  mean  value  consistent  with  the  volume  heating.  In  the  lower 
portion  of  the  burning  region,  the  temperature  increases  as  r  *  0.  The 
center-line  temperature  is  approximately  10  percent  higher  than  the 
boundary  temperature.  The  air  temperature  drops  rapidly  for  y  >  1  and 
falls  to  within  7  percent  of  the  ambient  at  about  300  m  above  the 
ground.  At  that  point,  a  well-developed  vertical  flow  is  established, 
representing  the  early  stages  of  the  convection  column  (see  Fig.  9). 

Figure  11  presents  contours  of  constant  perturbation  pressure. 

In  the  lower  portion  of  the  Cuming  region,  radial  pressure  gradients 
are  readily  observed.  The  maximum  pressure  drop  observed  in  the  field 
is  0.118  psi.  Above  the  burning  zone,  the  radial  pressure  gradients 
are  much  smaller  and  the  perturbation  pressure  decreases  rapidly  to 
zero . 

Figures  12  and  13  display  horizontal  and  vertical  velocity  pro¬ 
files,  respectively,  for  the  turning  region.  The  fire-wind  velocity 
peaks  (18.3  m/sec)  at  the  inlet  to  the  fire  zone.  The  induced  fire 
winds  decrease  rapidly  above  y  “  1.  As  y  ■+■  10,  the  induced  horizontal 
velocity  decreases  toward  zero,  indicating  that  lower  order  terms  in 
an  expansion  for  u  are  needed  to  describe  the  subsequent  entrainment 
of  ambient  air  into  the  column.  Classical  plume  theories  relate  that 
velocity  to  the  vertical  velocity  at  the  column  center  line.  As  in 
the  Flambeau  simulation,  profiles  of  vertical  velocity  (Fig.  13)  are 
"top  hat"  in  form.  The  maximum  velocity  predicted  at  the  top  of  the 
turning  region  is  4.19  m/sec. 
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Fig.  9--Velocity  field  for  Hamburg  firestorm 
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V.  DISCUSSION 


The  model  developed  here  describes  the  velocity  and  thermodynamic 
fields  generated  by  a  large  urban  fire.  The  analysis  focuses  on  the 
turning  region,  which  includes  the  burning  zone  and  the  region  below 
the  established  free-convection  column.  Such  an  approach  allows  esti¬ 
mates  to  be  made  of  the  conditions  necessary  for-  shelter  design  and 
of  the  environment  facing  survivors  and  rescue  workers.  In  addition, 
the  results  of  the  analysis  provide  insight  into  the  complex  inter¬ 
actions  occurring  in  Che  turning  region. 

The  analytical  procedure  employs  asymptotic  methods  to  simplify 
the  basic  equations  of  momentum  and  energy  conservation.  Such  simpli¬ 
fication  allows  one  to  consider  the  principal  force  balances  while 
neglecting  lower  order  effects.  A  f inite-voliune  heat  source  is  used 
to  model  the  combustion  processes,  and  large  changes  in  temperature 
and  density  are  allowed.  A  one-parameter  eddy-viscosity  model  is  used 
to  describe  the  turbulent  stresses,  and  a  graybody  approximation  em¬ 
ployed  to  model  hot  gas  and  smoke  radiation.  Jump  conditions  are 
derived  to  describe  the  rapid  changes  in  physical  quantities  at  the 
fire  periphery.  Those  conditions  effect  model  problem  closure,  allow¬ 
ing  the  induced  fire  winds  to  be  calculated  directly,  without  extensive 
far-fieid  computations. 

The  use  of  one- parameter  models  to  represent  the  combustion  pro¬ 
cesses,  turbulent  exchanges,  and  radiative  heat  transfer  simplifies 
the  analytical  formulation,  but  the  models  of  course  only  approximate 
the  complex  physical  phenomena.  Further  modeling  could  improve  the 
approximations,  but  the  current  understanding  of  the  turning-region 
hydrothermodynamics  may  not  support  a  more  rigorous  analysis.  In  our 
sample  simulations  of  the  Flambeau  fires  and  the  Hamburg  firestorm  of 
1943,  observational  data  are  used  whenever  possible  to  specify  the 
model  parameters.  Parameters  that  cannot  be  defined  from  available 
data  are  allowed  to  vary  over  a  range  of  values. 

A  closed-form  solution  to  the  relevant  boundary-value  problem  is 
developed  for  the  restricted  case  of  weak  heating.  That  solution 
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reveals  tha  structure  of  the  general  solution,  and  concisely  shows  how 
the  heating,  production  of  buoyancy,  pressure  gradients,  induced  fire 
winds,  and  formation  of  the  initial  free-convection  column  are  related. 
Those  relationships  are  also  evident  in  the  sample  numerical  simulations 
of  the  strongly  heated  Flambeau  and  Hamburg  fires.  At  the  fire  per¬ 
iphery,  the  temperature  increases  rapidly,  generating  steep  pressure 
gradients  and  high-velocity  fire  winds.  The  inrush  is  turned  up  by 
buoyancy  all  along  the  burning  zone,  and  forms  the  base  of  a  free- 
convection  column  in  the  upper  portion  of  the  turning  region.  The  gas 
temperature  increases  slightly  toward  the  center  of  the  fire,  but  falls 
rapidly  with  height  above  the  burning  zone. 

Finally,  the  model  provides  a  flexible  framework  useful  for  fur¬ 
ther  studies  of  the  large-fire  environment.  For  example,  the  solution 
defining  the  gas  dynamic  state  provides  information  prerequisite  to 
estimation  of  species  concentrations.  Leading-order  versions  of  the 
conservation  equations  for  each  species  can  be  developed  [Small, 

Larson,  and  Brode,  19811  and  solved  to  provide  predictions  of  large- 
fire  oxygen  and  noxious  gas  levels.  The  present  analysis  may  also  be 
extended  to  consider  time -dependent  heat  release  rates  as  long  as  they 
vary  slowly  over  time  periods  characteristic  of  the  resulting  convec¬ 
tion.  Part  II  presents  an  example  of  such  an  extension. 
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SUMMARY 


This  report  analyzes  the  velocity,  temperature,  and  pressure  dis¬ 
tributions  of  a  large  area  fire.  The  region  of  interest  is  the  burning 
area.  The  analysis  considers  the  effects  of  massive  heating  by  combus¬ 
tion;  the  results  presented  are  applicable  to  the  lower  surface  layer. 
Predictions  of  the  basic  street-level  environment  can  be  made  and  above¬ 
shelter  air  temperatures  can  be  defined.  Coupled  with  an  analysis  of 
debris  distributions  and  residual  nuclear  radiation  (fallout  or  Induced 
activity)  in  an  attacked  city,  the  velocity  and  temperature  predictions 
will  aid  in  defining  the  environment  facing  rescue  forces. 

This  analysis  differs  considerably  from  previous  treatments.  No 
attempt  is  made  to  extrapolate  results  valid  for  a  free-convection 
column.  The  theory  presented  directly  relates  the  heat  release  of  com¬ 
bustion  and  size  of  the  burning  area  to  calculated  velocities  and 
thermodynamic  properties.  The  induced  fire  winds  depend  on  the  size 
of  the  fire  and  the  combustion  rate  and  are,  to  leading  order,  inde¬ 
pendent  of  the  properties  of  the  convective  plume. 

The  results  are  restricted  to  large  fires  (small  fires  generate  a 
different  balance  of  forces  than  those  considered  here).  The  induced 
fire  winds  increase  in  intensity  with  heat  release  or  fire  size.  Com¬ 
pressibility  effects  limit  the  maximum  velocities  attained;  however, 
a  continuous  spectrum  of  fire-wind  velocities  is  predicted,  ranging 
from  very  low  levels  (corresponding  to  the  experimental  Flambeau  fires) 
to  high-velocity  winds  such  as  generated  by  the  1943  Hamburg  firestorm. 

It  is  clear  that  the  fire  is  the  forcing  agent,  and  external  in¬ 
fluences  are  not  needed  to  produce  the  devastating  velocities  encountered 
in  the  firestorms  of  World  War  II.  If  the  fire  persists  long,  swirl 
may  slightly  affect  the  high-velocity  winds  predicted  here. 

Sample  calculations  are  presented  for  the  largest  Flambeau  fires 
and  for  the  Hamburg  firestorm.  The  results  agree  with  the  documented 
observations  and  provide  insight  into  the  detailed  hydrothermodynamic 
interactions  occurring  in  a  large  area  fire. 
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A  nuclear  weapon  explosion  can  Ignite  fires  over  much  larger  areas 
than  those  covered  by  World  War  II  firestorms.  The  theory  developed  In 
this  report  has  been  applied  to  large  fires  such  as  may  occur  in  modern 
urban  areas  subject  to  megaton-yield  explosions.  Those  results  are 
presented  in  Part  II  of  this  report. 


